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Abstract 

We obtain the strong asymptotics of polynomials p„(A), A G C, orthogonal with respect to measures 
in the complex plane of the form 

where s is a positive integer, f is a complex parameter and dA stands for the area measure in the plane. 
Such problem has its origin from normal matrix models. We study the asymptotic behaviour of Pn(A) 
in the limit n, A —>■ oo in such a way that n/N —> T constant. Such asymptotic behaviour has two 
distinguished regimes according to the topology of the limiting support of the eigenvalue distribution 
of the normal matrix model. If 0 < \t\^ < T/s, the eigenvalue distribution support is a simply 
connected compact set of the complex plane, while for |tp > T/s the eigenvalue distribution support 
consists of s connected components. Correspondingly the support of the limiting zero distribution of 
the orthogonal polynomials consists of a closed contour contained in each connected component. 

Our asymptotic analysis is obtained by reducing the planar orthogonality conditions of the poly¬ 
nomials to an equivalent system of contour integral orthogonality conditions. The strong asymptotics 
for the orthogonal polynomials is obtained from the corresponding Riemann-Hilbert problem by the 
Deift-Zhou nonlinear steepest descent method. 
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1 Introduction 


We study the asymptotics of orthogonal polynomials with respect to a family of measures supported on the 
whole complex plane. To set up the notation for the general case, let Pn(^) denote the monic orthogonal 
polynomials of degree n such that 


Pn(A)pm(A)e = hn,N^r, 


n, m = 0,1, 2,... 


( 1 . 1 ) 


/c 


where W : C —>■ M is called the external potential, TV is a positive parameteij^ dA(A) is the area measure in 
the complex plane and hn^jv is the norming constant. The external potential is assumed to have sufficient 
growth at infinity so that the integrals in (1.1) are bounded. 

Planar orthogonal polynomials satisfying (1.1) appear naturally in the context of normal matrix models 


m where one studies probability distributions of the form 

1 




■'n^N 


— 




( 1 . 2 ) 


^For simplicity we use the simpler notation Pn(A) instead of pn,N{^), even though the orthogonal polynomials depend on 
the value of the scaling parameter N. 
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where A/at is the algebraic variety of n x n normal matrices 


Mn = {M-. =0}cMat„xn(C), 

and dM is the volume form induced on Afn which is invariant under conjugation by unitary matrices. 
Since normal matrices are diagonalizable by unitary transformations, the probability density (1.2) can be 
reduced to the form |43] 


1 


|Ai Xj 






i<j 


where Xj are the complex eigenvalues of the normal matrix M and the normalizing factor called 

partition function, is given by 


i<j 

The statistical quantities related to eigenvalues can be expressed in terms of the orthogonal polynomials 
Pn{X) defined in In particular, the average density of eigenvalues is 

j=o 

The partition function Zn^N can be written as a product of the normalizing constants 


Zn,N 


Y{hj^N- 

i=o 


While the asymptotic density of eigenvalues can be studied using an approach from potential theory [49] , 
the zero distribution of orthogonal polynomials remains an open issue for general potential weights despite 
general results in m- The density of eigenvalues Pn,N{X) converges (in the sense of measures) in the limit 


n —)■ oo , —)■ oo , 


N 

n 


1 

f ’ 


to the unique probability measure p* in the plane which minimizes the functional [251135] 


(1.3) 


■^(m) = j^loglA-ry| ^dp{X)dp{T]) +^ 


W{X)dp{X). 


(1.4) 


The functional I{p) in (1.4) is the Coulomb energy functional in two dimensions and the existence of a 
unique minimizer is a well-established fact under mild assumptions on the potential W(A) [49] . If W is 
twice continuously differentiable and its Laplacian AW is non-negative, the equilibrium measure is given 
by 

dp*{X) = AWiX)xDiX)dA{X), 


where xd is the characteristic function of the compact support set D = supp(/i*). Sub-leading order 
corrections to the behaviour of the eigenvalues distribution Pn,N{X) as n,iV —>■ oo and fluctuations in the 
bulk and at the boundary of the support D have been considered in [gisiiiaiiiis]. 
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The measure can be also uniquely characterized by the Euler-Lagrange conditions 


1T(A) 

T 


- 2 


J log |A - r]\dfi*{ri) + £d > 0 


(1.5) 


for all values A S C with equality in (1.5) on the support]^ of The Lagrange multiplier £d is called 
the (generalized) Robin constant. It is a non-trivial problem to determine the shape of the support set 
D. In some cases this problem is called Laplacian growth. When the potential W{X) is real analytic the 
boundary dD is a hnite union of analytic arcs with at most a finite number of singularities [48) . see also 

[ 33 ]. 

There is only a handful of potentials W (A) for which the polynomials Pn{^) can be explicitly computed. 
The simplest example is W{X) = |Ap for which the orthogonal polynomials Pn{X) are monomials of degree 
n, the constants hn^N and the average density of eigenvalues Pn,N can be computed explicitly in terms of 
the Gamma function. The matrix model associated to this potential is known as the Ginihre ensemble 
|8ILI30j and the density of eigenvalues Pn,N converges to the normalized area measure on the disk of radius 
Vt centered at the origin. In general, for radially symmetric potentials W = IT(|A|), the orthogonal 


polynomials are always monomials and in the limit (1.3), the eigenvalue distribution is supported either 
on a disk or an annulus by the Single-ring Theorem of Feinberg and Zee [28] , whose rigorous proof can be 
found in [52] ■ The correlation functions in this case have been studied in [T5] . 

The harmonic deformation of the Gaussian case W (A) = | Ap + (harmonic) has been intensively studied. 
In particular the potential W{X) = |Ap — t{X^ + A^) is associated to the Hermite polynomials for |t| < 1 /2. 


In the limit (1.3) the distribution of eigenvalues is the normalized area measure on an ellipse, while the 


distribution of the zeros of the orthogonal polynomials is given by the (rescaled) Wigner semicircle law 
with support between the two foci of the ellipse [23] . 

The normal matrix model with a general deformation W{X) = |Ap + Re(P(A)) where P{X) is a 
polynomial of a fixed degree has first been considered in the seminal papers m and [M] (see also the review 
article [56)1 . where the connection with the Hele-Shaw problem and integrable structure in conformal 
dynamic has been pointed out. More general potentials have been considered later in |55j . For such 
potentials, however, the matrix integrals have convergence issues in the complex plane and therefore a 
natural cut-off has been introduced in the work of Elbau and Felder [ 53 ] . In [ 21 ] the polynomials associated 
to such deformed Gaussian potentials have been studied and it was argued that the Gauchy transform of 
the limiting zero distribution of the orthogonal polynomials coincides with the Cauchy transform of the 
limiting eigenvalue distribution of the matrix model outside the support of the eigenvalues. Moreover, it 
was also conjectured in [21] that the zero distribution of the polynomials Pn(A) is supported on tree-like 
segments (the mother body, see definition below) inside the compact set (the droplet) that attracts the 
eigenvalues of the normal matrix model. 

For the external potential W{X) = |Ap -I- Re(tA^), Bleher and Kuijlaars [TT] defined polynomials 
orthogonal with respect to a system of unbounded contours on the complex plane, without any cut-off and 
which satisfy the same recurrence relation that is asymptotically valid for the orthogonal polynomials of 
Elbau and Felder. They then study the asymptotic distribution of the zeros of such polynomials confirming 
the predictions of [ 21 ]. Similar results were obtained for the more general external potential [ 55 ] W{X) = 
|Ap — Re(t|A|^), k > 2 and |t| sufficiently small so that the eigenvalue distribution of the matrix model 
has an analytic simply connected support. Cases in which the eigenvalue support has singularities were 
analyzed in [39] and [6]. In particular in the work [^ the external potential IT(A) = |Ap —2clog |A —a| with 
c and a positive constants, has been studied and the strong asymptotics of the corresponding orthogonal 


^To be precise, the equality on the support is valid only up to a set of capacity zero m- 
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polynomials has been derived both in the case in which the support of the eigenvalues distribution is 
simply connected (pre-critical case) or multiply connected (post-critical case) and critical transition was 
observed (see also m^)- We remark that in the work [6], differently form the previous works, the zeros 
of the orthogonal polynomials do not always accumulate on a curve that corresponds to the mother-body 
of the domain where the eigenvalues of the normal matrix models are distributed. 

In this work we study the strong asymptotic of the polynomials Pn(A) orthogonal with respect to a 
density where the external potential is of the form W{\) = jAp'* — tA'* — tA^ with A € C, s a 

positive integer and t G C \ {0}. By a simple rotation of the variable A the analysis can be reduced to the 
case of real and positive t. Therefore, without loss of generality, we may and do assume that t G K+, that 
is, 

W(A) = |Ap®-t(A®-bA") AgC, sgN, t>0, (1.6) 

and the associated orthogonality measure is 

e"^(l^l'""*(^'’+^'’))dA(A) A G C, s G N, t > 0. (1.7) 

Note that the potential W{X) has a discrete rotational Z^-symmetry. It was observed in [S] (see also [55]) 
that if a potential W (A) can be written in the form 

W(A) = ^Q(A*), 

s 

the equilibrium measure for W can be obtained from the equilibrium measure of Q by an unfolding 
procedure. In our particular case 

Q(u) = s|up — st(u + u) = s|u — — st^ 


corresponds to the Ginibre ensemble so that the equilibrium measure for the potential Q is the normalized 
area measure of the disk 

fr 

\u-t\ = 4, G = \ —, 


where T has been defined in (1.3). The equilibrium measure for W turns out the be equal to 

= ^|Ap(«-i)y^(A)dAl(A), 

where xd is the characteristic function of the support set 

D:={Xg C, |A® -t| < tc}. 


( 1 . 8 ) 


(1.9) 


We observe that for t < tc the equation ( |1.9[ ) describes a simply connected domain in the complex 
plane with uniformizing map from the exterior of the unit disk in the ^-plane to the exterior of D given 

by 


/(0 = «li + ^i 


with inverse 


F(A) = /-1(A)=G^A 
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For t > tc the domain defined by the equation (1.9) consists of s connected components which have 
a discrete rotational symmetry. For s = 2 the domain D is called Cassini oval. The boundary of D, 
namely dD can be identified with the real ovals of the Riemann surface S defined by S := {(X,ri) G 
C^, (A® — — t) = which has genus (s — 1)^. Such Riemann surface does not coincide with the 

Schottky double of I? for s > 2. 

The boundary of the domain D can be also expressed by the equation 


A = 5(A), 5(A) = t + 


ti 


A« -t 


( 1 . 10 ) 


The function 5(A) is analytic in a neighbourhood of dD and it is called the Schwarz function associated 
to dD (see e.g. UHl). 


Remark 1.1. The domain D is a quadrature domain with respect to the measure d/r*. Indeed for 
any function h{X) analytic in a neighbourhood of D one has, by applying Stokes’ theorem and the Residue 
Theorem, 

f h{X)dfi*{X) = ^ [ h{X)X^-^{S{X)ydX = Vcfch(Afc) 

Jd 2Tntj, Jqjj 


where 5(A) is the Schwarz function (1.10), Ck = - and Xk=t‘e 

- s 


1.1 Statement of results 

The goal of this manuscript is to determine pointwise asymptotics of the polynomials Pn{X) defined in 
orthogonal with respect to the weight ( |1.7[ ) in the two cases 

• pre-critical: t < t^', 

• post-critical: t > t^- 

The Z^-symmetry of the orthogonality measure ([TtI) is inherited by the corresponding orthogonal poly¬ 
nomials. Indeed the non-trivial orthogonality relations are 


f p„(A)A^*+'e-'^^(^)dA(A) =0, j = 0,..., fc - 1, 

JC 


where k and I are such that 


n = ks + I, 0 < I < s — 1, 

i.e., the n-th monic orthogonal polynomial satisfies the relation 

Pn{e^X) = e^pniX). 

It follows that there exists a monic polynomial of degree k such that 


p„(A) = yql!\xy. 


( 1 . 11 ) 
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Therefore the sequence of orthogonal polynomials {pn(A)}^o can be split into s subsequences labelled 
by the remainder I = n mod s, and the asymptotics along the different subsequences can be studied via 
the sequences of reduced polynomials 




Z = 0,1,..., s — 1. 


By a simple change of coordinates it is easy to see that the monic polynomials in the sequence 
are orthogonal with respect to the measure 




S — I — 1 , N 

7 := - e [0,1), 


( 1 . 12 ) 


(1.13) 


namely they satisfy the orthogonality relations 

f g['^(u)u^|ur2T'e-^(l“l'-‘“-‘“)d7l(u) =0, j = 0,...,fc-l. 

Jc 

As a result of this symmetry reduction, starting from the class of measures (1.71, it is sufficient to consider 
the orthogonal polynomials with respect to the family of measures (1.12). It is clear from the above 
relation that for I = s — 1 one has 7 = 0 and the polynomials ^\u) are monomials in the variable 
{u — t), that is, 

= {u-tf. 

It follows that the monic polynomials pks+s-i{^) have the form 

Pks+s-l{\)=\^-\\^ -tf. 


Remark 1.2. Observe that the weight in the orthogonality relation (1.13) can be written in the form 

|y|-27g-JV(|up-tu-ifi) _ log |u|)^ 

and it is similar to the weight with W(u) = |up — clog \u — a| with c > 0 studied in However 

in our case it turns out that c = —2j/N < 0, so the point interaction near a = 0 is repulsive and the 
asymptotic distribution of the zeros of the polynomials (1.13) turns out to be substantially different from 
the one in m- 

Dehne 

tl 


and for r > 0 the function 


ZO 


= log(t - y) + ^ -logrt + -— 

tzo Zq 


Let us consider the level curve 

Cj. := |a G C, Re (^r(A) =0, |A® — t| < rt 


}■ 


(1.14) 

(1.15) 


The level curves Cr consist of s closed contours contained in the set D, where D has been defined in (1.9). 
For these curves we consider the usual counter-clockwise orientation. Define the measure D associated 
with this family of curves given by 


diy = 


1 


27ris 


#r(A), 


(1.16) 


and supported on Cr- 
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Lemma 1.3. The a-priori complex measure di) in {1.16} is a probability measure on the contour Cr defined 
in (|1.15|) for 0 < r < fi. 

I 

Let us denote by v{pn) the zero counting measure associated with the polynomials Pm namely 

’^{Pn) = - V Sx, 
n ^' 

P„(A)=0 

where 6x is point distribution with total mass one at the point A. 


Theorem 1.4. The zeros of the polynomials Pn{X) defined in (1.1) behave as follows 

• for n = sk s — 1, let ui = . Then <«, wt =,..., = are zeros of the polynomials pks-i-s-i with 

multipicity k and X = 0 is a zero with multiplicity s — 1. 

• for n = ks -\-1, I = 0,..., s — 2 the polynomial Pn{X) has a zero in X = 0 with multiplicity I and the 
remaining zeros in the limit n, N ^ oo such that 


N = 


n — I 
T 


accumulates on the level curves Cr as in (1.15) with r = 1 for t < tc and r = for t > t^. Namely 
the curve C on which the zeros accumulate is given by 


C : 


(t - A®) exp ( ^ 


t 


pre-critical case 0 < t < tc 


post-critical case t > tc, 


(1.17) 


with I A® — t\ < zot. The measure v in (1.16) is the weak-star limit of the normalized zero counting 
measure Vn of the polynomials pn for n = sk -\-1, I = 0,..., s — 2. 












Figure 1.1: The blue contour is the boundary of the eigenvalue distribution support D defined in (1.9) 
relative to the normal matrix model and the red contour is C defined in (1.17) which describes the support 
of the limiting zero distribution of the polynomials p„(A). Here s = 3 and t < (left figure) and t > tc 
(right figure). 


Remark 1.5. We observe that the curve (1.17) in the rescaled variable z = 1 — takes the form 

1 


C : 


1-2 

j 

ze ^0 



zoe-‘ 


-L-l 


pre-critical case zq > 1, 
post-critical case 0 < zq < 1- 


(1.18) 


ti 


with zo = ^ and \z\ < zg. The curve C is similar to the Szego curve {z S C : |ze^ = 1, |z| < 1} 


that was first observed in relation to the zeros of the Taylor polynomials of the exponential function 
and coincides exactly with such curve in the critical case Zq = 1. The Szego curve also appeared in the 
asymptotic analysis of the generalized Laguerre polynomials, see e.g. The curve ( |1.18[ ) is 

the limiting curve for the zeros of the polynomials 


T^kiz) := ^^^-ql!\-t{z - 1 )) 


where q^l\u) have been defined in (1.11). 


From Theorem 1.4 and (1.5), the following identity follows immediately. 


Lemma 1.6. Given the potential W{X) in (1.6), the measure v in (1.16) with r = 1 fort < tc and r = 
fort > tc, and the contour C defined in (1.17), the equation 

di>{v) 




c ^ ^ ’ 


X^C, 


(1.19) 


defines the boundary of a domain which coincides with D defined in (1.9). The measure p,* in (1.8) of 
the eigenvalue distribution of the normal matrix model and the measure v of the zero distribution of the 
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orthogonal polynomials are related by 


f dp.*{T]) 

D A-77 



A € C\D. 


( 1 . 20 ) 


Remark 1.7. The identities (1.191 and (1.201 in Lemma 1.6 are expected to hold in general for a large 
class of normal matrix models. It has been verified for several other potentials (see, for example 


We also observe that for the orthogonal polynomials appearing in random matrices, in some cases, 
the asymptotic distribution of the zeros is supported on the so-called mother body or potential theoretic 
skeleton of the support D of the eigenvalue distribution. We recall that a measure v is a strong mother 
body for a domain D with respect to a measure p, if 


1 ) f log |A — r]\dp{r]) < f log |A — ri\dv{r]), for A S C with equality for A outside D; 


2) n > 0 and supp v C supp p; 

3) the support of v has zero area measure; 

4) the support of v does not disconnect any part of D from D‘^. 


If the measure v has only the property 1), 2) and 3) it is called weak mother body. The problem of 
constructing mother bodies is not always solvable and the solution is not always unique m- 

Concerning the explicit examples appearing in the random matrix literature, for the exponential weight 
W{X) = |Ap +Re(P(A)) where P{X) is polynomial, the support of the zero distribution of the orthogonal 
polynomials is a strong mother body of the domain that corresponds to the eigenvalue distribution of the 
matrix model (see, e.g., in [23 [M [13 In contrast, in the model studied in m and also in the 

present case, the support of the zero distribution of the orthogonal polynomials does not have property f) 
and therefore it is a weak mother body of the set D. 

The proof of Theorem |1.4| is obtained from the strong and uniform asymptotics of the polynomials 
Pn(X) in the whole complex plane which is obtained by characterising the orthogonal polynomials Pn{X) 
via a Riemann-Hilbert method. 

Uniform Asymptotics. In the next theorem we describe the strong and uniform asymptotic of the 
polynomials Pn(X) in the complex plane. We distinguish the pre-critical case t < tc and post-critical case 
t > tc. We first define the function 


0(A) = 


0r— jSq ('^) 0 <C t tc 

0r— l(A) t 


(pre-critical) 

(post-critical), 


( 1 . 21 ) 


where ((riX) is given by (1.141. 


Theorem 1.8 (Pre-critical case). For 0 < t < tc the polynomial Pn{X) with n = ks 1, I = 0,..., s — 2, 
7 = € (0,1), have the following asymptotic behaviour for when n, N ^ 00 in such a way that 

NT = n-l: 


(1) for X in compact subsets of the exterior of C one has for any integer M >2 


( 1 . 22 ) 
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(2) for X near C and away from A = 0, 

p-k(p{\) 


PniX) = 


1 + 


/ji+7 I r(-7) zo 


1\ ^ t 


x^ r A« 


^ /I 
+ 0 - 


(1.23) 


where ^(A) has been defined in (1.21) and r( 2 ;) = {z — 1)!; 

(3) for X in compact subsets of the interior of C and away from A = 0, 


Pn{X) = X‘ 


.-rk {-tf+^ 1 
[ r(-7) A« 


1 - 


1 

Zo 


-1-7 


o 


(4) for X in a neighbourhood of X = 0, we introduce the function w{X) = (j){X) + 27ri if A® € C_ and 
w{X) = <^(A) i/A® S C+. Then 




-k<j>{X) / /I 

(iI;(A))^-^^(vI/i2(fcu>(A)) + 0f- 


(1.24) 


where the {l,2)-entry of the matrix 'll is defined in (3.22). 


We observe that in compact subsets of the exterior of C there are no zeros of the polynomials Pn{X). 
The only possible zeros are located in A = 0 and in the region where the second term in parenthesis in 


the expression (1.23) is of order one. Since Re^(A) is negative inside C and positive outside C it follows 
that the possible zeros of Pn{X) lie inside C and are determined by the condition 


1 


t 


Ite0(A) = -(l + 7)t|i + llogl 


1 

-1-7 

t 

1 — 


1 - ^ 

^0 


A® 


|A®-t|<t. (1.25) 


The expansion (1.23) shows that the zeros of the polynomials Pn{X) are within a distance O (l/fc^) from 
the level curve (1.25). This level curve approaches C defined in ( 1.17|) at a rate O (logfc/fc). 
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Figure 1.2: The blue contour is the boundary of the support of the eigenvalue distribution defined in (1.91 
and the red contour C defined in (1.17), is the support of the limiting zero distribution of the orthogonal 
polynomials p„(A). The yellow contour is given by (1.25) in the pre-critical case and (1.27) in the post- 

critical case and lies within a distance of order O j from the contour C. The dots are the zeros of 
the polynomial Pn(A) for s = 3, n = 285, I = 0 and t < tc on the right and t > tc on the left, respectively. 

For t > tc the polynomials Pni^) with n = ks + I, I = 0,..., s — 2, have the following asymptotic 
behaviour. 

Theorem 1.9 (Post-critical case). For t > tc the polynomials Pn{^) with n = ks +1, / = 0,..., s — 2, and 


1 = 


S — l — 1 


€ (0,1) have the following behaviour when n, N ^ oo in such a way that NT = n — I 


(1) for A in compact subsets of the exterior of C one has 

p„(A) = A'(A* - + t{zo - 1))^ (^1 + O , 

t^ 

with zo = -I; 

F 

(2) for A in the region near C and away from the points A® = t(l — Zq) one has 

<((A® -t)A®)T' 


p„(A) = A'(A* - + t{zo - 1))^ 1 - 


g-fc0(A) ^^2 


C (A® -f t(zo - 1))2 t'+i 



(1.26) 


with 4>{X) defined in (1.21) and the constant c defined in (4.24); 

(3) for A in compact subsets of the interior region of C one has 


Pn{\) = a ' 


"0 t^zl ftzo''^ 


kh+'i 


A®7 


+t{zo-l)p+^ 


Ol I 
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(4) in the neighbourhood of each of the points that solve the equation A® = t(l — zq) one has 
p„(A) = A'(A- - ("<-A - + O (i 


where U{a;^) is the parabolic cylinder function (Chap. 13, Jl^) satisfying the equation = 

+ afjU and w^{X) = —(^(A) — 2Tri for A'* S C_ and w^{X) = —^(A) for A'* S C+. 


We observe that in compact subsets of the exterior of C the polynomials Pn{X) have zero at A = 0 with 
multiplicity 1. The other possible zeros are located in the region where the second term in parenthesis in the 


expression (1.26) is of order one. This happens in the region where Re 4>{X) < 0 and Re 4>{X) = O {\ogk/k), 
namely in a region inside the contour C and defined by the equation 


Re(^(A) = -(1 + 7)^^ 


t^zl |(A--f)A-r 

|c| |A" + t(zo-l)| 27 +i 


(1.27) 


The expansion (1.26) shows that the zeros of the polynomials Pn{X) are within a distance O (l/fc2)from 
the level curve (1.27). This level curve approaches C defined in ( 1.17| ) at a rate O (logfc/fc). 

The proofs of Theorem |1.8| and Theorem |1.9| are obtained by reducing the planar orthogonality relations 
of the polynomials Pn{X) to orthogonality relations with respect to a complex density on a contour. 
More precisely, the sequence of polynomials Pn(A) can be reduced to s families of polynomials ql^\x^), 
Z = 0,..., s — 1 as in (1.11) with n = ks + 1. The orthogonality relations of the polynomials 


are reduced to orthogonality relations on a contour. We then reformulate such orthogonality relation as 
a Riemann-Hilbert problem. We perform the asymptotic analysis of the polynomials ql!\x) using the 
nonlinear steepest descent/stationary phase method introduced by Deift and Zhou [55] and successfully 
applied in Hermitian random matrices and orthogonal polynomials on the line in the seminal papers, 
m .|20j. See |5T| for an introduction to this method, with a special emphasis on the applications to 
orthogonal polynomials and random matrices. It would be also interesting to explore the asymptotic for 
the polynomials p„(A) using the 9-problem introduced in |5S] . 

The zeros of Pn{X) accumulate along an open contour as shown in Figure [Hi The determination of 
this contour is a first step in the analysis. Riemann-Hilbert analysis in which a contour selection method 
was required, also appeared in the papers 1,110]. 

The Riemann-Hilbert method gives strong and uniform asymptotics of the polynomials PniX) in the 
whole complex plane. The asymptotic behaviour of the polynomials Pn{X) in the leading and sub-leading 
order can be expressed in terms of elementary functions in the pre-critical case t < t^, while in the post- 
critical case t > tc we have used, in some regions of the complex plane, parabolic cylinder functions as 


in mmM- The proof of Theorem |1.4| is then deduced from the strong asymptotic of the orthogonal 
polynomials. 

The paper is organized as follows: 

• In Sect.j^ by using the symmetry properties of the external potential, we reduce the orthogonality 
relations of the polynomials Pn{X) on a contour and recall how to formulate a Riemann-Hilbert 
problem for such orthogonal polynomials m- 

• In Sect.j^we obtain the strong and uniform asymptotics for the polynomials Pn{X) in the pre-critical 
case t <tc. 
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In Sect. 4 we find the strong and uniform asymptotics for the polynomials Pn{^) in the post-critical 


case t > tc and we prove Theorem 1.4 


In the critical case t = the set D in (1.9) that supports the eigenvalue distribution of the normal 


matrix models has a singularity in A = 0. The corresponding asymptotics of the orthogonal polynomials 
seems to be described in terms of Painleve IV equation as in m and this situation is quite different 
from the generic singularity that have been described in terms of the Painleve I equation |41).|39j. Such 
problem will be investigated in a subsequent publication. 

Furthermore, from the norming constants of the orthogonal polynomials, one may consider the problem 
of determining the asymptotic expansion of the partition function Zn^N in the spirit of m- Also this 
problem will be investigated in a subsequent publication. 
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2 The associated Riemann—Hilbert problem 


In this section we set up the Riemann-Hilbert problem to study the asymptotic behaviour of the orthogonal 
polynomials PnW- As seen above, the analysis of PnW can be reduced to that of the polynomials (A^^) 
introduced in (1.11). The polynomials q^\u) are characterized by the symmetry-reduced orthogonality 
relations 

=0, j = 0,...,fc-l. (2.1) 

Jc 


2.1 Reduction to contour integral orthogonality 


The crucial step in the present analysis is to replace the two-dimensional integral conditions (2.1) by an 


equivalent set of linear constraints in terms of contour integrals. In what follows it will be advantageous 
to perform the change of coordinate 

u = —t{z — 1), z G C, 
and characterize q^!'^ (u) in terms of the transformed polynomial 

T^kiz) := - 1 )) 

in the new variable z. The polynomial 7rfc(z) is also a monic polynomial of degree k and it can be 
characterized as follows. 

Theorem 2.1. The polynomial iTkiz) is characterized by the non-hermitian orthogonality relations 


j> '^k{z)z- 


e 


-Nt^z 


Z-1 


dz = 0 j = 0,1,... ,fc - 1 , 7 6(0,1), 


( 2 . 2 ) 
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where 'E is a simple positively oriented contour encircling z = 0 and z = 1 and the function j 
analytic in C\[0,1] and tends to one for \z\ —>■ oo. 



Proof. In order to prove the theorem we first show that the orthogo¬ 
nality relation for the polynomials q\.{u) on the plane can be reduced to 
an orthogonality relation on a contour. To this end we seek a function 
Xj{u,u) that solves the i9-problem 

(2.3) 


Figure 2.1: The contour S 


Having such a function, for any polynomial q{u) one has 
d[q{u)xj{u,u)du] = q{u)di[Xjiu,u)du A du 


where d denotes the operation of exterior differentiation. If such Xjiu,u) exists, one can use Stokes’ 
theorem and reduce the planar orthogonality relation to an orthogonality relation on a suitable contour. 
The equation (2.3) has a contour integral solution 


Xjiu,u) = 

Jo 


fgj-'r+i ^ 

d- 

-u) 

' 1 

1 1 



7 

]X[j-~t+i ' 

u) 

' C 

r(j-7 + i) 1 


t 

Xfj-z+i I 

U 


1 ^Ntv 


nNu{u — t) 


^dr 


^Ntu 


{u — 


r(j - 7 + 1 ) - 


r^ '^e ''dr 


/ Nu{u—t) 


7 ^Ntu 


uj (u - ty+^ . 


1 - O 


M —>■ oo 


It follows that for any polynomial q(u) the following integral identity holds: 

= 1 ™ / q{u)Xj{u,u)du 

2l i?->oo J\u\=R 

= - lim / qlu) \gAu) - O 

= — ® q{u)Gj{u)du 
2* 


where R and i?o are sufficiently large and Gj{u) = (l ~ u)'* (ix-t)j+i • follows that for any 

polynomial q{u) the following identity is satisfied: 


(z(u)u^|u|-2T'e-^(l“l'-‘“-‘“)dH(u) = 


7rr(j-7-H) I 
iV^-T'+i 2TTi 


^Ntu 

q(u) -r^TT 

\u-ty+^ 


1 - - 
u 


du 
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where 7 € (0,1), j is an arbitrary non-negative integer, and S is a positively oriented simple closed loop 
enclosing u = 0 and u = t. Making the change of coordinate u = —t{z — 1 ) one arrives to the statement 
of the theorem. Q.E.D. 


2.2 The Riemann—Hilbert problem 

Our aim is to study the behaviour of the polynomials 7r^,(z) in the limit fc —>■ cx) and —>■ 00 in such a 

way that for n = ks + I one has 

r>o. 


Let 


and introduce the function 


In terms of the weight function 


T 


T 

tc = — and zq = — 
s 


V{z) =\ogz-\ -. 

zo 


wu{z) := 


the orthogonality relations (2.2) can be written in the form 


j) nk{z)z^Wk{z) dz = 0 j = 0,1,... 


k-l. 


( 2 . 4 ) 


In the limit k —)■ cx), we distinguish two different cases: 

• pre-critical case zq > 1, corresponding to 0 < t < tc > 

• post-critical case zq < 1, corresponding to tc <t . 

Our goal now is to characterize the polynomial TTk{z) as a particular entry of the unique solution of a 
matrix-valued Riemann-Hilbert problem. Let us first define the complex moments 


j) z^Wk{z)dz 


where, for simplicity, the dependence on k is suppressed in the notation. Introduce the auxiliary polynomial 


nfc-i(z) 


1 

det [^i+j]o<ij<k-i 


det 


1^0 

Vl . . 

Vk-1 

1^1 

V2 ■ • 

Vk 


k’k -2 


k' 2 k -3 

1 

z 



( 2 . 5 ) 


Note that Ilfe-i is not necessarily monic and its degree may be less than k— 1 : its existence is guaranteed 
just by requiring that the determinant in the denominator does not vanish. 
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Proposition 2.2. The determinant det[^'i+j]o<i,j<fc-i does not vanish and therefore Wk-i{z) is well- 
defined. 

Proof. We have 


det[z/i+j]o<ij<fc-i = det 






z-1 


dz 




„-Nt^z 


0<i,j<k-l 
7 


Z-1 


dz 


0<ij<k-l 


where the last identity has been obtained by the reflection of the column index j —>■ k — 1 — j. Due to 
Theorem [2JJ we have 


[ Tr{z){z — iy\z — 1\ '^'^e dA{z)=f^ 

Jc 


-2,-27^r(j-7+l) 1 


iVi-7+1 2m 


7r(z)- 


^-Nt^z 


Z-J + 1 V 2 — 1 


dz , 


an hence the second determinant is given by 

, . . r-Nt^z 

det 



r 

C-i) H 

= det 


zyz-iy\z-l\-^'^e-^*"^^^"dA{z) 


Finally, the determinant on the right-hand side is strictly positive because 


k-l 

n 

-■ j =0 


2j^2i+27-2^i-7+l 


r(j-7 + i) 


( 2 . 6 ) 


det 


zyz-iy\z - ir2T'e-^*'l"l'd7l(z) 


= det 




> 0 , 


where the equality follows from the fact that the columns of the two matrices are related by a unimodular 
triangular matrix, while the inequality follows from the positivity of the measure. Finally, since r( 2 ;) has 
no zeros (and no poles since j — 7 + 1 > 0), the non-vanishing follows from (2.6). Q.E.D. 

Define the matrix 


7rfc(z) — 


1 f T^k(.z') 


Wk{z')dz' 


Y(z\- 2mJ^z'-z 

-27rillfc_i(2:) - / —;- Wk{z )dz 

Jt. z' - z 


(2.7) 


It is easy to verify that the matrix Y (z) is the unique solution of the following Riemann-Hilbert problem 
of Fokas-Its-Kitaev-type m- 


1 . Piecewise analyticity: 

Y{z) is analytic in C \ E and both limits Y±(z) exist along E, 


( 2 . 8 ) 


2. Jump on E: 


Y+{z) 



Wk{z)\ 

1 J 


z € E, 


(2.9) 
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3. Behaviour at z = oo: 


Y{z) ={l + 0 


, z -)> (X). 


The last relation is obtained by noticing that if Ilfc-i is given by (2.5) then the entry Y 22 


matrix T(z) in (2.7) satisfies 


nfe i(z') 


z — z 


Wk{z')dz' = z 1 + 0 -)), z—>■ 00 . 


2.3 Initial undressing step 

In order to simplify the subsequent analysis, we define the following modified matrix: 


f(z) := Y(z) 1 - 


1\ 2' 


zec\(su[0,i]). 

This matrix-valued function Y (z) satisfies the following Riemann-Hilbert problem: 
1. Piecewise analyticity: 


Y (z) is analytic in C \ (S U [0,1]) 


2. Jumps on E and [0,1): 


f+(z) = y_(z) 


2 ^-kV(z) 

0 1 


z G E, 

ze (0,1), 


with V(z) as (2.4). 

3. Behaviour at z = 00 : 


Y{z)=[l + 0[\ 


Z^’^^ , Z 00 . 


4. Endpoint behaviour at z = 0 and at z = 1: 

Y{z)z-^'^^ =0{l) z^O, f(z)(z-l)T‘"^ =0(1) z^l. 

The polynomials (z) is recovered from Y (z) as 

/ 1 \ 2 

T^k{z) = Yii(z) 1 - 


( 2 . 10 ) 

(z) of the 


( 2 . 11 ) 
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3 Asymptotic analysis in the pre-critical case 

In order to analyse the large k behaviour of Y[z) we use the Deift-Zhou nonlinear steepest descent method 
|22j . The first step to study the large k behaviour of the matrix function Y(z) is to make a transformation 
Y{z) —)■ U{z) so that the Riemann-Hilbert problem for U{z) is normalised to the identity as \z\ —)■ oo. 
For this purpose we introduce a contour C homotopically equivalent to E in C \ [0,1] and a function g(z) 
analytic off C. Both the contour C and the function g{z) will be determined later. We assume that the 
function g(z) is of the form 

giz) = J log{z - s)dv{s), (3.1) 

where di'{s) is a positive measure with support on C such that 


With this assumption clearly one has 


dv{s) = 1. 


g{z) = log z + O {^z az |z| —>• 00 , 
where the logarithm is branched on the positive real axis. 


(3.2) 


3.1 First transformation F —)■ f/ 

Since C is homotopically equivalent to E in C \ [0,1] we can deform the contour E appearing in the 
Riemann-Hilbert problem for Y io C. Define the modified matrix 

U{z) = e-^-(^/2)-3F(^)e-fc9(D-3gfeW2)<T3 ^ g c \ (C U [0,1]), 

where £ is a real number, to be determined below. Then U(z) solves the following RHP problem 

1. Piecewise analyticity: 

U{z) is analytic in C \ (C U [0,1]). 

2. Jump discontinuity on C: 

k{g+-g-) „k{g^+g--l-V)\ 

0 j ■ 


U^{z) = U.{z) 


3. Jump discontinuity on (0,1): 


U+{z) = , 


4. Endpoint behaviour at z = 0 and z = 1: 

I7(z)z-i^^ =0(1) z^O, U(z)(z-l)i^^ =0(1) z^l. 
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5. Large z boundary behaviour: 


t/(z) =/ + 0 

The polynomial 7rfc(z) is determined from U{z) by 

7rfe(z) = ( 1 

3.1.1 The choice of g-function 


(3.4) 


(3.5) 


In order to determine the function g and the contour C we impose that the jump matrix (3.3) becomes 
purely oscillatory for large k. This is accomplished if the following conditions are satisfied: 


g+{z)+g.{z)-e-V{z)=0 
Re(g+(z) - 5 -(z)) = 0 


z eC. 


(3.6) 


Next we show that we can find a function g and a contour C that satisfy the conditions (3.6). For the 
purpose we use the following elementary result in the theory of boundary value problems. 

Lemma 3.1 f[29). p. 78). Let L be a simple closed contour dividing the complex plane in two regions 
and D_ where = Int(L) and D_ = Ext(L). Suppose that a function '0(C) defined on L can be 
represented in the form 

0(C) = 0+(C) + 0-(C), CeL 

where ip+iC) is the boundary value of a function '0+(z) analytic for z € D+ and '0-(C) is the boundary 
value of a function ijj-{z) analytic for z S D_ and such that 0_(oo) = 0. Then the Cauchy integral 




^(0 

c-^ 


dC, 


can be represented in the form 


$+(z) =-0+(z), forzeD+, 
<i)_(z) = —'if’-{z), for z G D-. 


The boundary values of the function <i> on the two sides of the contour L then satisfy 

CGL. 


4>+(C) + $-(C) = 0+(C)-V'-(C) 
4>+(C)-$-(C) = 0+(C) + V'-(C) 


Now we apply Lemma 3.1 to the function g'{z) that satisfies the differentiated boundary condition 

5+(C) + 5-(C) = f Tj CgC. 
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The functions ip±{z) of Lemma b.ll are therefore the function g'iz) 

\ \ zo z 

takes the form 


= Sri 


1 _ i_ 

^ - ^-^dT = 

Z — T 


1 

Zo 

1 

z 


z € Int(C) 
z € Ext(C). 


(3.7) 


so that the measure diy in (3.11 is given by 


dvlz) = -— f-1 dz, z €C. 

2-ki \z Zo/ 


Integrating the relation (3.7) and using (3.2) one has 


9{z) = 


-z £ Int(C) 

Zo 

logz z G Ext(C) 


(3.8) 


(3.9) 


where I is an integration constant and logz is analytic in C\K+. Performing the integral in (3.1) for a 
specific value of z G C, say z = 0 and deforming C to a circle of radius r one can determine the value of t. 

r 

I = logr-, r > 0. 

Zo 


The total integral of dv in (3.8) is normalized to one on any closed contour containing the point z = 0. 
However we have to define the contour C so that dv{z) is a real and positive measure along C. To this end 
we introduce the function 


(priz) = 


—2g{z) + V{z) + logr —— = logz —^— logr + — z G Int(C)\[0,r] 

Zo Zo Zo 

2g{z) — V{z) — logr + — = logz —^— logr + — z G Ext(C)\(r,oo) 
Zo Zo Zq 


(3.10) 


Observe that (pr{z) is analytic across C, namely (pr+(z) = (pr-{z), z G C and that 

, / N (pr+{z) + (pr-{z) 

(pr{z) = --- =-g+{z) + g-{z), zGC. 


(3.11) 


The next identity follows in a trivial way 


dv{z) = —d<pr{z). 

2tti 


Imposing the relation (3.6) on (3.11) one has 

Re((;ir.(z)) = log |z| - _ log r + — = 0. 

Zo Zq 

This equation defines a family of contours which are closed for |z| < r (easy to verify). Since the function 
r 

— logr H-for r > 0 has a single minimum at r = zq and diverge to +oo for r —>■ 0 and r —>■ oo, it is 

Zo 

sufficient to consider only the values 0 < r < zq. 
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We define the contour Cr associated to r as (see Figure 3.1) 


Cr = {z S C : Re((/)r( 2 :)) = 0, \z\ < Zq}, 0 < r < zq- 


(3.12) 



Figure 3.1: The family of contours Cr for 
three values of r: r = zq, r = 1 and r < 1. 
For r = Zq, the region \z\ > zq is plotted 
as well. 


\x\ J- 

Since dx^Q4>r\y=o = —y-> 0, for |a:| < zq, Re(l)r{z) is 

X Zq 

negative inside Cr- Note that the point z = 0 lies inside Cr 
since Re (j){z) —>■ —oo for \z\ —>■ 0. Furthermore Cr intersects 
the real line in the points z = r and 0 S (—r, 0). Indeed 
Re (j){—r) > 0 which shows that the point z = —r lies outside 
Cr- 

Lemma 3.2. The a-priori complex measure v in 
probability measure on the contour Cr defined in 
0 < r < Zq . 

Proof. By the Residue Theorem the measure diz has total 
mass 1 on any contour Cr- Since Ret)),. = 0 on the contour 

Cr it follows that the measure dv = - d(f>r is real on the 

2Tri 

contour Cr- In order to show that the measure is positive on 
the contour Cr we introduce the variable 


(3.8) is a 
(ra for 


fir 


= e^ 


On the contour Cr we have that \fir\ = 1, so that in the ^^.-plane the contour Cr is mapped to a circle of 
radius one and the map fir = e'^'' is a univalent conformal map from the interior of Cr to the interior of a 
circle (52]. We have 

1 dfir 
2TTi fir ’ 

which shows that the measure dv in the variable fir is a uniform measure on the circle and therefore dv 
is a positive measure on each contour Cr- Q.E.D. 

We are now ready to prove Lemma |1.3| and Lemma |I.6| announced in the Introduction. 


Proof of Lemma 1.3 By using the Residue Theorem, it is straightforward to check that dh is a 
probability measure on any contour Cr defined in (1.15). In the post-critical case such curve Cr has s 
connected components Cfi j = 0,..., s — I, each encircling exactly one root of the equation A'* = t- In the 
pre-critical case we denote with the same symbol C-j. the connected components of Cr \ {0} together with 
{ 0 }. 

With the appropriate choice of the branch of the root, any can be parametrised by the map 


A(z) := (-t(z- I)); 
implies 


which respects the orientations of C and C^ . Since 4>{z) = ^(A(z)), equation (1.16) 
di'{X{z)) = -dh'{z) - 


This relation, together with Lemma 3.2 where we proved that dv^z) is real and positive on Cr, implies 
that also dv is real and positive on any C^ and hence on the whole C-r- Q.E.D. 
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Proof of Lemma 1.6 By using the Residue Theorem we first calculate the r.h.s. of (1.19) obtaining 

., 8-1 


1 


V 


1 


f dv(ri) 

Id ^ ~ V 27rf 7(5 A — 7? 77 ® — t ' ZQt J 


+ 


dr] = 


y 


-t' 


(3.13) 


Therefore (1.19) is equivalent to 

sX^-\y-t) A'*-! 


T 


or (A®-t)(A®-t) = - 


which is equivalent to the equation for the boundary of D defined in (|1.9[). Next we show that the l.h.s. 


obtains 


of (1.20) is equal to its r.h.s. By using Stokes’ Theorem, the relation ( 1.10[ ) and the Residue Theorem one 

?7®“^77®d77 1 f r]‘^~^S{r])dr] A'*”^ 

A® - t 

Q.E.D. 


dfi*{r]) 


1 


In 27r5t2 


'dD 


X — rj 


1 

2-Kit}. 


c JdD 


X-r] 


which coincides with (3.13). 


3.1.2 Choice of the contour 


^ 

- 

1 

1 

1 

1 

t 0*- 


\ • s 



^ ^0 \ 

^-- 

^r<l 


C 

\ 

\ 

1 

J 

/ 


We now argue that the relevant contour on which the zeros 
of the orthogonal polynomials Kk{z) accumulate in the pre- 
critical case is given by the level r = 1. The family of contours 
0 < r < 1 can be immediately ignored because in this case E 
has to be deformed to two contours E ~ Cr<i U C as shown 
in Figure 3.2 On the contour C we have Re((/)) > 0 and 


Figure 3.2: The contour Cr<iUC that is ho¬ 
motopic to E. The region where Re < 
0 is coloured in red. 


it is not possible to perform any contour deformation to get 
exponentially small terms in the jump matrix ( |3.3| ) as A: —> 00 . 

For 1 < r < zo the original contour E can be deformed 
homotopically to the contours Cr- In this case the asymptotic 
expansion of the matrix entry Yii(z) = Kk{z) as k ^ 00 , does 
not give any sub-leading contribution. For this reason we 
discard this case and we omit the corresponding asymptotic 
analysis. The only possible case remains r = 1. In this case 
the analysis as fc —>■ 00 to the matrix entry Fii(z) = Kk{z) 


which will performed in subsections 3.2 to |3.5[ gives leading 
and sub-leading terms. The comparison of these two terms 
enable us to locate the zeros of the polynomial Kk{z) on a contour that lies within a distance O (logfc/fc) 
from the contour C^^i- 

So for the reasons explained above, we are going to perform the asymptotic analysis of the RH-problem 
(3.3|-(3.4| by deforming the contour E to the contour Cr—i- For simplicity we denote this contour by C: 


C = ( z e C : log |z| — -I- — = 0, |z| < 1 

zo zo 
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1 


Also, for simplicity, the function (j)r=i will be referred to as (j) below: 

(3.14) 


= \ogz- — + —, 
Zo Zq 


where logz is analytic in C\(—oo,0]. With this choice of contour, the 
jump matrices of the RH problem (3.3)-(3.4) for the matrix U{z) are 
summarized in Figure [373| 


Figure 3.3: The jump matrices 
for U{z) 


3.2 The second transformation U ^ T 

Consider two extra loops Ci and Cg as shown in Figure [T4l These define 
new domains Oq, fli, 112 and Iloo- Define the new matrix-valued function 


Z G Doo U Dn 


Tiz) = 


U{z) 

t/(^) Z&V12. 

Then T{z) satisfies the following Riemann-Hilbert problem: 


(3.15) 


1. Piecewise analyticity: 

T{z) is analytic in C \ (Cg U C U Ci U [0,1]) 

2. Jump discontinuities on St = Cg U C U Ci U [0,1]: 

T+iz) = T_{z)vt{z), Z&Y.T, 


where 


vt{z) 



g— 77ri(T3 


Z & C 

Z eCe 

Z G Ci 
ZG (0,1). 


3. Endpoint behaviour at z = 0 and z = 1: 

T{z)z-i'^^ =0{1) z^O, T(z)(z-1)^‘"=‘ =0(1) z^l. 


4. Large z boundary behaviour: 

T(z) = 1 + 0 , z —)■ 00 . 

Note that is analytic in C \ {0} . The important feature of this Riemann-Hilbert problem is that 

the jumps are either constant or they tend to the identity matrix as A: —>■ oo at an exponential rate. 
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Proposition 3.3. There exists a constant cq > 0 so that 


vt{z) = I + O {e as fc —>■ oo 


uniformly for z € Ce U Ci\Ui, where U\ is a small neighbourhood of z = 1. 


The proof of this proposition follows immediately from the fact that Re ( 5 ^( 2 ) < 0 for z S Int(C)\f//i and 
Re^(z) > 0 for z G Ext(C)\Z^i. Here Int(C) is the interior region bounded by C and Ext(C) is the exterior 
region bounded by C. 

From the above proposition it follows that 



3.3 The outer parametrix for large 2 ; 


We need to find a matrix-valued function P°°{z) ana¬ 
lytic in C \ (C U [0,1]) such that it has jump discontinu¬ 
ities given by the matrix v°°(z) in (3.16), namely 


Figure 3.4: The jump matrices for T(z) with 
contour Ey = Ce U C U U (0,1). 


P“(z) =Pf°(z) 





Z G C 

z e (0,1) , 


with endpoint behaviour 

P“(z)z-^'"^ = 0(1) z^O, P“(z)(z-1)5'"^ = 0(1) z^l, 


and large z boundary behaviour 


Define 

where 


P°°{z) = I + 0 [ - ) , z^oo. 


P^iz) ■.= P^{z)x-\z) 


X{z) := 


0 1 
-1 0 
/ 


z G Int(C) 
z G Ext(C) . 


(3.17) 


The A:-independent matrix P°°{z) has no jump on C and it satisfies the following Riemann-Hilbert 
problem: 
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1. Piecewise analyticity: 

2. Jump discontinuity on (0,1): 

3. Boundary behaviour at z = oo: 


P°° is holomorphic in C \ (0,1) 


P^{z) = I + 0\ - 


oo 


The solution to this Riemann-Hilbert problem is given by 


P°°{z) = 1- - 

' 2 


which leads to the particular solution 


P“(z) = 1 - 


1\ 


X{z) = 


1 - - 
z 

1-i 


0 1 

-1 0 


z € Int(C) 
z € Ext(C) 


(3.18) 


3.4 The local parametrix at 2 : = 1 

The aim of this section is to construct a local parametrix P°(z) in a small neighbourhood Ui oi z = 1 
having the same jump property as T for z near 1 and matching the outer parametric P°°{z) in the limit 
/c —)■ 00 and z € dlAi. Then the Riemann-Hilbert problem for P^{z) is given by 

P+(2) =P°(z)ut(z), zG^JiHEt, 


and 

P'’(z) = P°°{z){I + 0 ( 1 )) as fc —)■ 00 and z G dUi. (3.19) 

In order to build such local parametrix near the point z = 1 we first construct a new matrix function 
P(z) from P°(z): 

P(z)=P°(z)x(z)-iQ(z) (3.20) 


where 

( /l pk(p{z)\ 

(Q 2 ) 2 G Oq n zJi 

I z G ZJi \ Uq 

The matrix P(z) satisfies the following jump relations in a neighbourhood of 1: 


Q{z) := 


P+(z) = P_(z) 


,77ri (g7'^* _ g-7'^i^gfc0(z) 


0 


77r2 


z G (— 00 ,1). 


(3.21) 


Next we construct a solution to the so called model problem, namely a Riemann-Hilbert problem that 
has the same jumps as B. 
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3.4.1 Model problem 

Consider the model problem for the 2x2 matrix function 'I' ) analytic in C\K“ with boundary behaviour 


^+(C) 


'i'-(c) ( 

( I + O 





3 — 7712 




Defining 

we obtain the following Riemann-Hilbert problem for ^(^): 


^+{0 = ^-{0 (J 


(1 - e-2T'^*)(^T')+e«^ 
^ ^ 00 , 


C G (-00,0) 
^ 00 . 


^ G (-00,0) 


where the function is analytic in C\(— 00 ,0] and (^"'')+ denotes the boundary from the right with respect 
to the contour (—oo,0) oriented in the positive direction. This is an abelian Riemann-Hilbert problem, 
and therefore easily solvable by the Sokhotski-Plemelj formula: 


/ 1 - 6 - 27 " r {C)+e<dC\ 

'f'(?)=| 2^* A- C-e 1 


(3.22) 


so that 

In particular 

g-277ri _ 2 ^ 

4-1 = 


4/(0 = 


LXJ T- 

i=i ^ 


A = 


0 1 
0 0 


as ^ 00 . 


27ri 


JR- 27rz Jo 


r(-7) ■ 


(3.23) 


(3.24) 


3.4.2 Construction of the parametrix 


We are now ready to specify the parametric P^{z) which follows from combining (3.231 with (3.20), that 
is 

P°{z) = Eiz)^{kw{z))Qiz)-^xiz ), 


where Q(z) and x(z) are defined in (3.21)) and (3.17) respectively and E{z) is an analytic matrix in a 
neighbourhood of Ui and w{z) is a conformal mapping from a neighbourhood of 1 to a neighbourhood of 

0 . 

The conformal map w(z) is specified by 


f (j){z) + 2Tri z GUi r\ C_ 

1 cj){z) z G n C_|_ . 
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We observe that 


w{z) =(^l-iyz-l)-^{z-lf + 0 {{z - 1)3) z ^ 1 . 

The matrix E[z) is obtained from condition (|3.19 1 which, when combined with (3.24) gives 


P°°{z){P°{z))-^ =P°°{z)x{z)-^Q{z)'i!{kw{z))-^E{z)-^ = 

1 A 


= P^{z)x{z)-\kw{z))-i^yi-y^^^^^PO{k-yE-\z), fc^oo, zeaWi 


(3.25) 


where we use the fact that Q{z) I exponentially fast as fc —>■ oo, and z G dUi. From the above 
expression it turns out that the matrix E{z) takes the form 


E{z) = P°°{z)x{z)~\kw(z))-i'^^ = 1 - - {kw{z))-i'^E 


(3.26) 


We observe that the function E{z) is single valued in a neighbourhood ofZ^i. Indeed the jumps of (z— 1 ) 3'^3 
and w(z)“ 5‘^3 cancel each other. 

From the expression (3.25) and (3.26) the matching between P°(z) and P°°{z) takes the form 

^M-l 






(3.27) 


1=1 


= /- 1 - 


' M-l 




as fc —>■ 00 


for z G dUi, where 7 G (0,1). 


3.4.3 Riemann-Hilbert problem for the error matrix R 


We now define the error matrix R in two regions of the 
plane, using our approximations to the matrix T. Set 


R{z) 


T{z){P^^\z)) \ 
T{z) (P“(z))-\ 


Z G 7/1 , 

everywhere else. 


(3.28) 


The matrix R is piecewise analytic in C with a jump 
across the contour C_r = Ci U Ce U dUi given in Fig. |3.5| 


RH problem for R 

1 . P is analytic in C \ Cij, 



Figure 3.5: The jump contour structure Cn 
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2. For z S Cr, we have 
with 


vr{z) 


R+{z) = R^{z)vr{z), 

P^{z)vt{z){P^{z))-^ zeCnXdUi 

P“(z) zedUi. 


(3.29) 


3. As z —>■ oo, we have 


P(z) =/ + 0(z-i) . 


The jump matrices across the contour Cr \ dUi are all exponentially close to I for large k because vr 
converges exponentially fast to v°° defined in (3.16) and the product P“(z)u°°(z) (P“(z)) = I with 

P°°{z) defined in (3.18). The only jump that is not exponentially small is the one on dU\ since from 
(3.27) for any integer M > 2 we have 

u«(z) = P-(z)(P°(z))-i 

1 


= /- 1 - 


'M-1 


{kw{z)) ^ ^ 




^ {kw{z)y 


0{k 


-M\ 


with the shift matrix A defined in (3.23). By employing the notation 

1\^ 




z j {w{z)y~^'f ’ 




we can rewrite the matrix vr{z) in the form 


1 


1 = 1 


ki 


+ 0(fc-“) A 


where, in particular, (3.24) implies 




1 - 


1 


(u>(z)) 


- 7-1 


By a standard perturbation theory argument one has the expansion 

1 


P(z) — J + 






which gives, using (3.291, (3.30) and (3.32) 


R^^\z)=R^y\z)+v^^\z)K. 


Therefore 


P(i)(z) = 



[ 4'\0dC' 

2TTi 7 

- 1 

1 


A 


A, 


(3.30) 


(3.31) 


(3.32) 
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where we observe that the function (z) has a simple pole in z = 1 with expansion 

- 7-1 


« ^ ’ ~ r(-7) z - 1 

By a simple residue calculation we obtain 


1 - 


R^^\z) = < 


(z- l)r(-7) 
1 


1 

1 - 

Zo 


1 - 

Zo 


- 7-1 


- 7-1 


Zo 


0 1 
0 0 


(l + 0(z-l)) . 


zGC\Ui 


0 1 
0 




(3.33) 


(z-i)r(-7) 

In general, given the structure of the jump matrix (3.301 the error matrix R^^'^{z) in (3.32) is of the form 


i?«w=(o y, 

namely, only the (l,2)-entry of the matrices R^^'^{z) is non-zero. From (3.28) and (3.32) one has 

+ E 

P°°(z) -h 


(3.34) 


T{z) = 


1 


' M-l 




y 


+ O {k-^) \ P^°\z) zGUi 


E + O ) P°°{z) everywhere else. 


(3.35) 


3 = 1 


where, in particular, the first term R^^'^{z) is given in (3.33). 


3.5 Proof of Theorem 1.8 asymptotics for J9„(A) for Zq> 1 

In order to obtain the asymptotic expansion of the polynomials Pn(A) for n —> oo, and NT = n — I we 
first derive the asymptotic expansion of the reduced polynomials Trk(z) for fc —> oo with n = sk + 1. For 
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the purpose we use (3.5), (3.15), (3.33) and (3.35) to obtain 


7rfc(z) = 1 -- 


[Uk{z\^ 
[Tk{z\^ 


= ( 1 - i ) < 


1 0 

gk(l>{z) ^ 


= e>^9{^) ( 1 - i ) < 


Tk{z) 


Tk{z) 

[i?(z)P“{z)]ii 

P(z)P“(^) f ® 


z G u r^o 

Z G 


{ 1 

oy 


ijj 


Z ^ ^^2 


11 


11 


{ 1 

oy 

y_g-k<p{z) 

ijj 


[R{z)P%z)\^^ 


f 1 

oy 

y^ktj>(z) 

ijj 


P(Z)P'^(Z) ( J 


z € r^oo u (iio \ ^i) 

z G \ IA\ 

Z € ^^2 \ ^1 

z € ^1 

z € r^i n z^i 

Z ^ ^^2 ^1 ■ 


From the above relation we obtain the expansions in the following regions. 

The exterior region Qoo- 


In this region, from (|3.34|) and (|3.35|), for any integer M > 2 we have 

1 


'^k 


(^)=gfes(^) _ ly +0 


^M+7 


= .Mi-i)7i + o(pP 


on any compact subset of floo- Therefore there are no zeros accumulating in this region. 

The interior region flo \ • 


nk{z)=e^s^^^ 


k'^+'f {z - l)r(-7) V Zo 


1 - 


- 7-1 


o 


fc2+7 


(3.36) 
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The leading term of the above expansion is of order O (k ^ ^), so there are no zeros in this region. 
The interesting region fli \ Ui. 

l^-7 / , X -7-1 

+ 0 


'^k 


(z) = _ 1 j 


^k4>{z) _^ (l z) ' (l- — 

*1+7 {z - l)r(- 7 ) 'v ^0 


1 


*2+7 


where is uniformly bounded on fli C {Re((/>) < 0}. 

The other interesting region 0 € 172 \*/i: 

„-k<p{z 


'^k 


(z) = (^1 - 


1 - 


* 7+1 


(z-l)r(- 7 )V Zo, 




where e is uniformly bounded on 172 C {Re((/)) > 0}. 

The region Ui. for z ^UiC] Ext(C) 

1 ~ u;(z)7 ( (1 “ 


® (zu;(z))7 


-7-1 


(3.37) 


(3.38) 


w{z)'^ — 


*7 


^'l 2 (*U’(z)) - 




for z G n Int(C) 

Mz) = 


(z-l)7 


(zw(z))7 


kMz) t N 7 «'l 2 (*w(z)) U;(z) 

e ’w[z) ^ — 


i'-i) 


- 7-1 


*7 


*1+7 ( 2 -i)r(- 7 ) 


+ ^Ri^) + 0{i 



where ^12 is the (1, 2)-entry of the matrix ^ defined in (3.22 1 and u^^(z) has been defined in (3.31). Here 
Int(C) and Ext(C) is the interior and exterior of C respectively. 

In order to obtain the asymptotic behaviour of the polynomials Pn{^) we use the substitution 


Pn(A) = (-t)''A‘ 7 rfe ( 1 - y ) , n = ks + l, 7 = 


s — 1 — I 


which gives the relations (1.22)-(1.24) in Theorem 1.8 Q.E.D. 

Proposition 3.4. The support of the counting measure of the zeros of the polynomials TTk{z) outside an 


arbitrary small diskUi surrounding the point z = 1 tends uniformly to the curve C defined in (1.18). The 
zeros are within a distance o(l/*) from the curve defined by 


Re(j){z) = -( 1 + 7 )^°^ + ^log 


/ 1 i+r 

1 1 

l |r(-7)||z-i|7+i 

Zo 


-7-1N 


(3.39) 


where the function (l>{z) has been defined in (3.14). The curves in (3.39) approachC at the rate O {log k / k) 
and lies in Int(C). The normalised counting measure of the zeros of TTk{z) converges to the probability 
measure v defined in (3.8). 
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Proof. Observing the asymptotic expansion (3.36) of 7rfe(z) in i}ao\Ui it is clear that 7rfe(z) does not have 
any zeros in that region, since z = 0 and z = 1 do not belong to floo \ ^i- The same reasoning applies to 
the region flo \Ui where there are no zeros of Trk(z) for k sufficiently large. 

From the relations (3.37) and (3.38) one has that in fli U O 2 using the explicit expression of g{z) 


defined in (3.9) 


TTkiz) = (1- - 


1 - 


o-k<l>{z) 




1-1 

Zo 


_(z- i)r(-7) 

The zeros of 7rfc(z) may only lie asymptotically where the expression 

p-k<j>{z) / 1 \ - 7-1 


-7-1 




1 - 


(z-i)r(-7) 


1 - 


z € 111 U ^21 


is equal to zero. Since 112 C {Re(0) > 0} and Hi C {Re((/)) < 0}, it follows that the zeros of 7rfe(z) may lie 
only in the region Hi and such that Re0(z) = O (logfc/fc). Namely the zeros of the polynomials 7rfe(z) lie 
on the curve given by ( |3.39| ) with an error of order O Such curves converge to the curve C defined 

(1.18) at a rate O (logfc/fc) (see Fig. 3.6). 



Figure 3.6: The zeros of 7r/j(z) for s = 3, Z = 0, t < tc, and k = 65, 80, 95. The red contour is C while the 
yellow contour is the curve (3.39). 


We now know that the zeros of 7rfe(z) accumulate on the curve C. We still need to determine the 
asymptotic zero distribution. Due to the strong asymptotic of TTk{z) we have from (3.36) 


lim ilog(7rfe(z)) = / log(z - 
k —^00 /c J C 


-C)dv{C)^ 


uniformly on compact subsets of the exterior of C. Furthermore C is the boundary of its polynomial convex 
hull. Then it follows that the measure dv in ( 3.8[ ) is the weak-star limit of the zeros distribution of the 

49j Chapter 3). 

Q.E.D. 


polynomials 7rfc(z) (see [44], Theorem 2.3 and 
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4 Post-critical case 

In this section we assume that 


0 < zo < 1 ■ 

To perform the analysis similar to the pre-critical case, one has to choose the appropriate contour from 
the homotopy class of the contour S that encircles the branch cut [0,1]. Since the point Zq is lying in 
(0,1), the family of curves , 0 < r < zg defined in (|3.12 ) are loops encircling z = 0 and crossing the 


real line in z = r. The regions where Ke(j)r{z) < 0 with 4>riz) defined in (3.10), are depicted in the figures 


for two values of the parameter r. 

) 

1 

1 

1 

1 

/ 

^ -^ 

1 

1 

1 

1 

1 

1 

/ 

/ 

{ ® * 

•zo • 

-N 

\ 

\ 

\ 

\ 

1 

1 

C ' 

s^r 1 

i 

1 

1 

1 

1 

1 

1 


Figure 4.1: 
Re(j)r{z) < 


The contour Cr for r = zg on the left and for 0 < r < zg on the right. 
0, with (j)r{z) defined in (3.10), is red. 


The region where 


In order to perform the asymptotic analysis of the Riemann- 
Hilbert problem (2.8), (2.9) and (2.10), we need to deform E to 


a homotopic contour U Co in such a way that the Re(j}r{z) is 
negative on Co¬ 


lt is clear from Figure 4.1 that the only possibility is to deform 
the contour E to the contour C^g UCo with Co as depicted in Fig. 4.2 
For simplifying the notation we define the following. 

Definition 4.1. Let 


- - - + 1 zgC\(0,+oo) , (4.1) 

Zq 


:= Kiz) = log 

4>+iz) — 4>-{z) = —27ri z G (0, +oo) . 


with 

and 



C :=C, 


with (('zo(z) defined in (3.10) andCzg defined in (3.12). 


According to (3.11) on the contour C we have 


Figure 4.2: The contour UCq. The 
region Re ft^g {z) < 0 is plotted in red. 


9 +{z)+g-{z) = V{z)+e 
9-e{z) - g-{z) = -(t>{z) , 
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where 


I = log(zo) - 1 ■ 


4.1 First transformation F i—)■ S' 


Since C U Co is homotopically equivalent to S in C \ [0,1] when zq < 1, we can deform the contour S 
appearing in the Riemann-Hilbert problem (2.8)-(2.10) for the matrix V(z) to CUCq. Define the modified 
matrix 

S(z) = 2 e C \ (C U Co U [0,1]), 

where V(z) has been defined in ( 2.11| ). Then the matrix S(z) is the unique solution of the following 
Riemann-Hilbert problem with standard large z behaviour at z = oo: 


1. S(z) is analytic for z G C \ (C U Co U [0,1]),. 

2. Jump discontinuities (with i^(z) as in ( |4.1[ )): 
ak4>(z)tT3 Z G C 

z GC„ 


S+{z) = S_{z) 


1 gfe<i>(z) 

0 1 

g-7^*<T3 ^ g (0, 1) 



3. Endpoint behaviour at z = 0 and z = 1: 

S{z)z-i'^^ =0{1) z^O 
S'(z)(z-= 0(1) z^l. 


Figure 4.3: The jump matrices for S'(z). The function 
(j){z) is defined in (4.1). 


4. Large z boundary behaviour: 


S(z) = /-|-0(-) , z—>- 00 . 


The orthogonal polynomials TTk{z) are recovered from the matrix S'(z) using the relation 

1 \ ^ 


^fc(z) = (^ 1 -^y Sn(^). 


4.2 The second transformation S T: opening of the lenses 

Consider two extra loops Ci and Cg as shown in Figure [T^ These define new domains Oq, Hi, 122 and 
12oo- Define 

z G 12 qq U 12q U 123 


T(z) = 


S{z) 

/ 1 

o\ 

S(z) 

^_gfe<A(2) 

1) 

S(z) 

f 1 

o\ 

g-fc0(z) 

1) 


Z G 122 . 
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Then this matrix-valued function has the following jump discontinuities: 


T+{z) = T_{z)vt{z), z 
0 1 


vt(z) = < 


-1 0 

1 0 

g-k<p(z) ^ 

1 0 

gfe0(z) 2^ 

1 gkcpiz) 

0 1 

-j-rrias 


Z G C 

Z G Ce 

Z G Ci 

Z G Co 
z G (0,1), 


where Yjt is the contour defined in the Figure 4.4 



Figure 4.4: The jump matrices for T[z) and the contour St = C U Ci U Ce U Co U (0,1). 
Endpoint behaviour: 

T{z)z-i'^^ = O {!) z^O, T{z){z - = O {!) z^l. 

Large z boundary behaviour: 

r(z) =1 + 0 , z —)■ 00 . 

Proposition 4.2. There exists a constant cq > 0 so that 

vt{z) = 1 + 0 as fc —> 00 

uniformly for z € Co U Ci U CffUza, where Uzq is a small neighbourhood of zq. 
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The proof of the proposition follows in a straightforward way by observing that by contruction (see 
Figure 4.2) Re^(z) > 0 for z G Ce\l^zo and Re(/>(z) < 0 for z G Ci UCo\Uzg. Therefore we have 

vt{z) -g v°°{z) as fc —>■ 00 


exponentially fast, where 


0 1 

-1 0 


as z G 


v°°{z) = < 


g-^*7<T3 as z e (0,1), 


as z G Ci U Ce U Co- 


The polynomials 7r/j(z) can be expressed in terms of T(z) in the following way: 

1 


^fc(z) = e'=9W(l--) 5n(z) = 


Tiiiz) 


= I 1 “ “ I 


= ( 1 - - 
z 


z G rioo u rio u ria 

z G rii 

Z G ^2 


Tii(z) z G rioo u rio u ^^3 

Tii(z) + e'=‘^(")Ti2(z) zGC!i 
Tn(z)-e-'='^(")Ti 2 (z) z G ^2 . 




T(z) 

^ J 


f 1 0 

T(z) 

_^-kcj>{z) 2 ^ 


(4.2) 


(4.3) 


4.3 The outer parametrix for large z 

Ignoring the exponentially small jumps and a small neighbourhood Uzg of zq where the uniform exponential 
decay does not remain valid, we are led to the following RH problem for : 

1. is holomorphic in C \ (C U [0,1]), 

2. satisfies the following jump conditions on C and (0,1): 

p|°°)(z) = p 1 °°^(z) 

[ z G (0,1) 

3. P^°°\z) has the following behavior as z —>■ oo, 

P^{z)=I + 0 


oo . 


0 

-1 


z G C 
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The above Riemann-Hilbert problem can be solved explicitly in terms of the piecewise defined matrix 
function 

V.):=| (-1 J) 

[ I z & Ext(C) . 

Define 

P“(z) := P^{z)x-\z) . 

The matrix P°°{z) has no jump on C and it satisfies the following RHP: 

1. P°° is holomorphic in C \ [0,1], 

2. Jump across (0,1): 


P“(z) =P“(z) 


2 g (0^ Zo) 

g-7^*<T3 2 g 


3. Large z boundary behaviour: 


P^{z) = I + 0 


z —>■ oo . 


A particular solution is given by 

P“(z) = 


which leads to the particular solution 


P“(z) 


'(z - Zq)^' 



{z - ZqY 

i<^3 

[z(z-l)J 


{z - Zo)^ 


'{z - Zof 

1^3 




z € Int(C) 
z € Ext(C) . 


(4.4) 


4.4 The local parametrix at z = zq 

The aim of this section is to construct a local parametrix P^{z) in a small neighbourhood U^q of Zq having 
the same jump properties as T(z) for z near zq and matching the outer parametric P°°{z) in the limit 
A: —)■ oo for z € dUzg ■ 


4.4.1 RH problem for P^{z) 

1. P'’(z) is analytic for z € Uz^ \ 

2. P)((z) = P^{z)vt{z) for z & Uzof^ St, 

3. for z G dUzo, we have 

P”(z) = P°°{z){I + o(l)) as fc —)• oo and z € dUzo- (4.5) 
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In order to build such local parametrix near the point z = zq we first construct a new matrix function 
B{z) from P°(z). Let us first define 


A(z) 


I Im(z) < 0 

e-T'"'"" Im(z) > 0, 


and the matrix Q as follows 

Qiz) = < 



for z e ria n 

for z e rio n u^a 

for z e (fii\fio) 
elsewhere. 


(4.6) 


(4.7) 


Then the matrix B{z) is defined from P'^{z) by the relation 

B{z) = P°(z)Q(z)e'=^(")"3/2A(z)-i. (4.8) 

The matrix B{z) satisfies the jump relations specified in Figure [T^ in a neighbourhood of Zq. In the next 
section we construct the solution of the so called model problem, namely a 2 x 2 matrix that has the 
same jumps as the matrix B. 


4.4.2 Model problem 


Consider the model problem for the 2x2 matrix function dt (^) 
analytic in C\{K U iK} with jumps and boundary behaviour 


vI/+(0 = 4/_(0u^(0 



as ^ oo , 



(4.9) 

- fT3 

(4.10) 


with dti 4'2 and 'I'a constant matrices (independent from ^) 
and where the matrix w$(C) is defined as 


u,t(^) = < 


1 1 - e-27^* 
0 1 
1 0 
g27'7ri 2^ 

g277rii ^ _ g -27 

0 e-27^* 

1 0 
-1 1 


c G M+ 

C G *+ 
e G K- 
c G -*+ 
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20 


Vo 1 ) 


5 


-2 



5-1 


Figure 4.5: The jumps for the matrix B 


The solution of the Riemann-Hilbert problem 
Recall the parabolic cylinder equation 


(4.9) and (4.10) 


is obtained in the following way, [22l [37l l46] . 




/ 
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/ = 0 , 


(4.11) 











(see, e.g., Chapter 19, [T]). This equation has a non-trivial solution U{a,^) associated to any a S C 
specified by the asymptotic behaviour 

+0 ^ ^ ^ oo, |arg(0| < |, (4.12) 

which is an entire analytic function of Three other solutions can be obtained by symmetry: 

U[a,-^), U{-a,iCj, U{-a,-i^). 


The relations among the above four solutions are 

r(i +a) 


U{-a, ±*C) = ±0 + tO) 

V 27r ' 


U{a,±0 = (e-“(“+5)/2^(-a,±iO + e“(“+5)/2z^(-a,T*C)) 

\/27r ^ / 


Moreover, solutions for a and a + 1 are connected by 


—U{a, + ^K{a, 0 + (® + 2 ^^^^ + 0 — 0 - 


(4.13) 

(4.14) 

(4.15) 


By using (4.13 )-(4.15), the solution of the Riemann-Hilbert problem (4.91 and (4.10) takes the form 




2-i<T3 


for ^ e 5±i 


'f'(C) = 


where 


eT‘2^f^(7- l,=piy2,c) ^ 

^ e±-^Z^(-7-|,-y2e) 

/3i, = -e— ^ 


(4.16) 


forCe5±2 


r(l-7)2T' 2/321 ■ 

Furthermore, from (4.12) one obtains the extra terms of the asymptotic expansion of 'I'(^) as ^ oo 

1 


(4.17) 




namely 


0 
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0 J e 


1 


-7(7-1) 0 


0 7(7 + 1)J 4^' 


1 




0 

(7-1)(7-2) 


(7+1)(7+2) 

421 




0 

^ 0 


^'2 = V 


-7(7-1) 0 

0 7(7 + 1) 


0 


(7+1)(7+2) 
^21 

(7-1)(7-2) q 


(4.18) 


(4.19) 
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4.4.3 Construction of the local parametrix 


From (4.8 1 and the previous section we are now ready to specify the form of the local parametrix at z = zq, 
by 

P°(z) = E(z)'^(^w(z))A(z)e-^^(^^^^Q(z)-^ 


where E{z) is an analytic matrix in a neighbourhood of the matrix dt has been defined in (4.16) 
and A(z) and Q(z) have been defined in (4.6) and (4.7) respectively. The function w{z) is a conformal 
mapping from a neighbourhood of zq to a neighbourhood of 0 and it is specified by 

w^iz) ;= 


—ipiz) — Sttz z G Uzo n c_ 

-^(z) zeitzo(^C+ 


We observe that 


w{z) = -4^{z - zo) - (z - zof + 0((z- Zof) 


\/ 2 z() 


Zo ■ 


(4.20) 


The matrix E{z) is obtained from condition ( |4.5[ ) which, when combined with (4.18), gives 
P“(z)(P°(z))-i = P“(z)Q(z)ei^(")'^^A(z)-i4'(Vfcn;(z))-^£;(z)-i 

+ O E~^{z), k ^ oo, z G dUz, 


= P“(z)x(z)-iA(z)-i(v^n;{z))-'^"^ (l-j 


0 

0 


'/kw{z) 


(4.21) 


where we used the fact that Q(z) —>■ x(z) ^ exponentially fast as /c —>■ oo, and z G dlAzg- From the above 
expression and (4.5), it turns out that the matrix E(z) takes the form 


E{z) = P°°{z)x{z)-^A{z)-\Vkw{z))-'^'^^ = 


{z - ZqY 
Lz(z-l)J 


A{z)-\'/kw{z))-^’^^ 


(4.22) 


We observe that the function E{z) is single valued in a neighbourhood otUzo- Indeed the boundary values 
of (z — zo)''"^^ and cancel each other. The boundary value of (z(z — !))(]( = (z(z — l))le^’^*'>', so 

that A(z)(z(z — 1 )) 5'^3 remains single valued in a neighbourhood of Zq . 

From the expression (4.21) and (4.22) the matching between P°(z) and P°°{z) takes the form 


P“(z)(P«(z))-i=P(z) J-i 


0 ^ 

P21 


2 I./3I2 J Vkw{z) 


+ O [k P ^{z), k^oo, z G dUzo, 


where 7 G (0,1). It is clear from the above expression that the (2, l)-entry of the matrix above is not 
small as fc —>■ 00 . For this reason we need to introduce an improvement of the parametrix. 


4.5 Improvement of the local parametrix 

In order to have a uniformly small error for A: 00 we have to modify the parametrices as in HU], [Sj, 

m- 

P“(z) := f/+ —P“(z) , 

\ Z — ZqJ 
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where C is a nilpotent matrix to be determined and 


P\z) := E{z) ( J 

\ Vkw(z) J 


where the matrix 't'l has been defined in (4.10) and 

m = (i+ 


c 


, E{z). 

Z- ZqJ 

With those improved definitions of the parametrices, P°°{z) and P°(z) have the same Riemann-Hilbert 
jump discontinuities as before but they might have poles at z = zq. Note that E{z) also has a pole at 
z = zq. However we can choose C in such a way that P^{z) is bounded in Zq. This is accomplished by 

-1 


C = -E{zo) 


where 


E{zo) = 


0 

^ 1,21 

Vkw'{zo) 

2zo 


E'izo) 


0 

'&l ,21 

Vkw'(zo) 


1 


0 ^ 


+ E{zo) I 'i’^ 2 iui''izo) , 

\ 2 \/k{w'k(zo)y 


^A:(l - Zq) 

From the above relation the matrix C takes the form 

0 0 
0 


2 ■ T 4zn — 1 

gTri^o-a E'{zQ)=^a^- -—- zE^Zq). 


6zo(l - zqY 


c = 


and 


1-^0 

2 zo 


e--*^zo(vI/i)2i25 = - 


r(i - 7 ) ( 1-^0 

'/YYt V Zq 


Zq. 


The improved parametrix gives the following matching between P°°{z) and P^{z) as k 

Zq £ dUzo '■ 


= E{z) ( I - ^1 - ^2 , 4 ' 24 'i + 4 'i 4'2 - 4'? - 4'3 


= Mz) 


I- 


Vkw{z) kw{zy 

0 (4'i)i2 

0 0 


k2w{zY 


1 


+ 0{k-^) ) ( 

'/kw{z) 


(4.23) 

(4.24) 
> 00 and 

E{z)-^ 


Vkw{z) 


(4'2)ii 0 

0 ( 4 ' 2)22 — 2 

kw{zY 


0 —(4'3)i2 

(4'i)2i(4'2)i1 — (4'3)21 0 

kiwiyzY 


0{k- 


E{z)- 


(4.25) 


which shows that P°°{z){P^{z)) ^ = I + 0 (l/fc“), a > 0, as fc —?> 00 when zq G dUzg for 7 G [0,1). 

4.5.1 Riemann-Hilbert problem for the error matrix R 
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We now define the error matrix R in two regions of 
the plane, using our approximations to the matrix 
T. Set 


R{z) = 


T(z) (P~(z))”\ z€C\U, 


'^Zq’ 

(4.26) 

The matrix R is piecewise analytic in C with a 
jump across (see Fig. 4.6 1 . 

Riemann Hilbert problem for R 



1. i? is analytic in C \ Sij, 

2. For z G E_r, we have 

with 


Figure 4.6: The contour U Ce U Cq U dU^g 

where Ci,Ce and Co are defined only in <C\Uzg- 


pc 


VRiz) = 


R+iz) = R_{z)vr{z), 


^{z)vTiz)(^P^{z)) ' ZGJ:R\dUz. 


(4.27) 


P' 


3. As z —>■ oo, we have 


(z) (p(°)(z))” 
R{z) = I + O (z-^) . 


Z G dUz, 


The jump matrices across the contour YiR\dlAzg are all exponentially close to I for large k because 
vt{z) converges exponentially fast to v°°{z) dehned in (4.2) for z G C\Uzo and the product 


pc 


’(z)^~(z) (p-(z)) 


-1 


= / 


c 


z - zo 


P“(z)n~(z) (P“(z)) 


-1 


C 


Z- Zo 


= I 


with P°°{z) defined in (4.4) and C in (4.23). The only jump that is not exponentially small is the one on 
dUzg- Indeed one has from (4.25) 


Ufl(z)=P“(z)(pO(z))-i = / + 


7e(z)2 


/ ckz 2 


2A;2+T'ri;(z)/32i 


-1 


\ 


z - Zo 
.2^27-1 


c‘/c‘' " ckz 2 

V(z-Zo)2 Z-Zo/ 


+ 


7(1 - 7) 

4kw{z)^ 


1 


fc 2 '''r(;(z)3 


O 

7(1 - 7) 

,4,di2e(z)2 


0{k-^z-j 

+ 0{k-^) 


O 


2ckz-h 
(z - Zo) 


0 

-1 


where 

e(z) = Pii(z)/c'^/^ 


(4.28) 


and we have substituted the explicit expressions of the matrix Ti, 4/2 and d'a as 
We have the following two cases depending on the value of 7 G (0,1). 


given by (4.19). 
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(a) 0 < 7 < i 


where 


(b) i < 7 < 1 


vr{z) = P“{z)(pO(z))-i = I + ^ + O (fc-i) , 


Vr{z) = - 


ie{zY 


2w{z)P21 V ® ^ 


0 1 


(4.29) 

(4.30) 


Vr{z) = P^(z){P^{z))-^ = I + 




(z) , V^^\z) 


n(3) 




+ 


{z) 


where 




4”(-) = - 


w{z)e{zY^ 


c(l - 7 ) 


fcs+'T 

( 1 - 7 ) 


O 


(k-’-i) . 


(4.31) 


0 0 

/92 i (z-zo)^ ' 2ri;(z)2(z-zo) ' 2/3i2w{z)^e{zy J \l 0 


7 c 


e(z)^ 


_, 7 ( 1 - 7 ) 

2^21 w(z)(z - Zo) 4 zi;(z )2 
7 e(z)^ /O 1 ' 


(^3 


(4.32) 


2w{z)f32i V® 


By the standard theory of small norm Riemann-Hilbert problems one has a similar expansion for R{z) in 
the large k limit. 


Case (a): 0 < 7 < 1 


r(i) 


(4.33) 


Compatibility of (4.27), (4.29) and (4.33) and the jump condition i?_|_ = R-Vr on dU^o gives the following 
relations. 

r!'1\z) = r!'1\z) + v^j^\z), z G dUzo 


with u^^(z) defined in (4.30). In addition R^^^z) is analytic in C\dUzo and ^^^^(oo) = 0. The unique 
function that satisfies those conditions is given by 




7e(^o)^ 


0 1 


z — zo 2 w'(zo)P 21 ® 


z G C \ lAz, 


= 


'^^r\z) 


leizoY 


Z — Zo 2w'( zq) (321 7 ^ 


0 1 


where the integral is taken along dUz^- Using the definition of P 21 , w'{z) and e(z) given in (4.17), (4.20) 
and (4.28) respectively one has 
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= < 


1 7^0 ^0 1 
0 0 


Z — Zq C 


v‘i\z) + 


1 7^0 fo 1 


Z — Zq C 


0 0 


Z £ C \ ly/z. 


z € lAz. 


with the constant c defined in (4.24). 
Case (b): | < 7 < 1 


R[z) =1+ ^ + ^ + ^ + O 




k 2 


{k-'-i) 


(4.34) 


Compatibility of (4.27), (4.31) and (4.34) and the jump condition i?+ = R-Vr on dUzo gives the following 
relations. 

R^+ {z) = (z) + (z), z G dUzo 

with v^^{z) defined in (4.32). In addition R^'^\z) is analytic in C\dUza and i?(*^(oo) = 0. 

Since has a third order pole at z = zq, the unique function that satisfies those conditions is 

given by 


R^^\z) = - 


ResA=zo v^R W R.esA=2o (A) ResA=zo (A - ^ (A) /O 0 


* 0 


Z G C \ k/zo ■ 


Z- Zq (z - Zo )2 (z - Zo )3 

Given the structure of the matrix P^{z), the matrix i?*^^^(z) does not give any relevant contribution to 
the orthogonal polynomials 7rk(z). 

Regarding R^‘^'> (z) one has 


_ ResA=zo(A-zo)r;)^^^(A) 
- - (z-zo )2 


,( 2 )/ 


z - Zo 

1 cje{zop /2e'(zo) w"(zo) 


z - Zo 2/32iry'(zo) V ^(^ 0 ) 2w'(zo) z - zq 

1 7(1 - 7 ) w"(zo) 1 A 

' as z G <L\Uz, 


0-3 


z-Zo4w'(zo)2 \ri;'(zo) z - Zq 


so that, using (4.17), (4.20), (4.24) and (4.28) one obtains 


R( 2 )(z) = < 


72:0 / 7(4 zo - 1) 2-7 3-7 

z - Zo V 3 z:o(l - Zo) 3 zo 2 (z-Zo) 


0 - 3 , 


( 2 ). X _ 7Z^0 f 7(4zo - 1) 2-7 3-7 A 

. ^ z - Zo V 3 zo(l - Zo) 3 zo 2 (z-Zo)y ^ 

In a similar way for R^^^ (z) we obtain 

^ _ResA=zoA(2!W _ 1 ie{zop 


z G C \ Llzi 


z G Liz 


z - Zo 


z - Zo 2r(;'(zo)/32i ® 


0 1 


Z G C \ Uzo , 
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so that 


1 l4 fO 1 


R^^\z) = < 


Z — Zq C 


v‘i\z) + 


0 0 


1 izi fO 1 
0 0 


Z — Zq C 


z G C\lAz. 


z € lAz. 


with the constant c defined in (4.24). 


4.5.2 Proof of Theorem 1.9; asymptotics for Pn{^) for 0 < zo < 1 

In order to obtain the asymptotic expansion of the polynomials Pni^) for n —)■ oo, NT = n- 
0 < zo < 1, we first derive the asymptotic expansion of the reduced polynomials 7rfc(z) as k oo. 


Using the relation (4.3) and the relation (4.26) one obtains 


7rfe(z)=e'=®(^) (1--) < 


i?(z)P“(z) 


11 


i?(z)P~(z) ( J 


P(z)P“(z)(_^i^(,) J 


P(z)P°(z) 


P(z)P‘'(z) ( J 


p(z)p‘'(z) J 


z € (Uoo U Uq U Ua) \ Lizo 
z € fli \ lAzo 
z G U 2 \ 

z G (Uqo u Uq u O3) n Uzq 
z € n 

z G ^^2 n hizQ ■ 


j 11 


The region (fioo U \ Uzo 
^fc(z) = 

The region Uq \ U^a 




= pfes(U ( 

z 




ki+^ I c (z-z,p+i+^[k))- 


with c defined in (4.24). 


I and 


(4.35) 
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The region Sli \ Uza 

■nk{z) = e 


= ^fes(z) { 


Oi I 


^24 


Z — Zq C 


(z - Zq)^ 

z(z-l) 


21 -7 


k 2 


+ < 


/ 7(4 zo -1) , 2-7 ^ 8-7 

Z-Zo A:V3zo(l-^o) 3zo 2{z - Zq) 


O 


1 




0<7< i 

H7<1 


with c defined in (4.24) and where we observe that Re(/)(z) < 0 in Oi. In a similar way we can obtain the 
expansion in the region 0,2 \ Uzg ■ 


The region O 2 \ U^q 

TTkiz) = e 


= pfeg(z) i ^0 
z 


X 7 / 

(M 

0 

Kf 

1 

1_ 

/ i fcH 7 C 

[z{z-l)\ 


Z- Zo 


OI i 


+ 


7^0 1 / 7(4 zo - 1) 2-7 

z - Zo k [3zo(l - zq) 3zo 
where now we observe that Rec/)(z) > 0 in 112■ 

The region l^/^g 


3-7 

2(z - Zo) 


a 


kl-y 


0<7< i 


< 7 < 1 


Using the relations (4.35), (4.16), (4.13) and (4.14) one obtains 

= + .eu,nExi(C) 

where U is the parabolic cylinder function that solves equation (4.11). For z S Uzg (7 Int(C) we have 


In order to prove Theorem |1.8| it is sufficient to use the above expansions, make the change of coordi¬ 
nates A® = —t{z — 1) and use for n = fcs -I- Z the relation 

Pn(A) = (-t)'"A'7rfc(l - A®/t). 

Q.E.D. 

Proposition 4.3. The support of the counting measure of the zeros of the polynomials TTk{z) for 0 < t < tc 
outside an arbitrary small disk Uzg surrounding the point z = zo tends uniformly to the curve C defined in 
(1.18). The zeros are within a distance o{\/k) from the eurve defined by 


Re(l){z) = - ( i -h7 


log(fc) 1 

k k 


log 



{z - Zq)'^ 


^ - 1 

1 


(4.36) 
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where the function (j){z) has been defined in ( |3.14 |. Such curves tends to C at a rate Oilogk/k). The 
normalised counting measure of the zeros of Trk{z) converges to the probability measure v defined in (3.8). 

Proof. As in the proof of Proposition |3.4[ it is clear from the above expansions of the polynomials (z) 
that there are no zeros in the region (floo U fls U Oq) \ sufficiently large k. Then we observe that 

the asymptotic expansion of the polynomials TTk{z) in the regions U ^2 \^ 2 o takes the form 

/ p-k(t>{z) 


7rfc(z) = z^ 


z- Zq 




1 - 


7^0 \{z:-zo)^ 


fcs+T' 


z{z — 1) J z — Zq 


-7 


+ 


7^0 1 f 7(4zo - 1) 

Z- Zok V3zo(l - Zq) 


2—7 3—7 


3zo 2(z - Zq) 

SO that we conclude that the zeros of 7rk(z) occur in the region where 


O 


kl-y 


0<7< i 


i <7< 1 


1 - 


g-fc<A(2) ^^2 


fc5+7 


(z - Zofi 


1 


z - Zq 


Iz(z-l)] 

is equal to zero. 

Since fl 2 C {Re((^) > 0} and C {Re((/)) < 0}, it follows that the zeros of 7rk(z) may lie only in 
the region Oi and such that Re(/>(z) = O (logfc/fc). Namely the zeros of the polynomials 7rfe(z) lie on the 
curve given by (4.36) with an error of order O (l/fc^). Such curves converges to the curve C defined (1.18) 
at a rate O (logfc/fc) (see Figure [4?7| ). 

k=65 k=80 ,, 



Figure 4.7: The zeros of TTk{z) ioi s = 3, I = 0, t > tc and k = 65, 80, 95. The red contour is C while the 
yellow contour is the curve (4.36). 


The proof of the remaining points of Proposition 4.3 follows the lines of the proof of Proposition 3.4 

Q.E.D. 


5 Proof of Theorem 1.4 


From the comments after the statement of Theorem 1.8 and Theorem 1.9 it is clear from (1.25) and (1.27) 


that the zeros of the polynomials Pn(A) accumulate in the limit n = ks + I ^ oo, with s and I given and 
Z = 0,..., s — 2, along the curve C defined in (1.17). In order to show that the measure k is the weak-star 
limit of the zero density we will show that 

- Odvif), 


lim -logp„(A)= / log(A-0d 

n—)-oo n 
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for A in compact subsets of the exterior of C (namely the unbounded component of C\C). Indeed using 
the relation 

Pn(A) = (-t)''A'7rfe(z(A)) 

one has 

lim -logp„(A)= lim ^ log(A^(-t)*^7rfc(2:(A))) = - log(-t) + - lim ^ log7rfc(z(A)) . 

n—>-oo n k—^oo kS + / 5 5 k—^oo k 

Using Proposition 3.4 and the relation diy{z{X)) = sdi){X) we get 

- log(-t) + - lim I log7rfc(2:(A)) = - log(-t) + - / log(z - 

5 5 k^oc k 5 5 J(2 


= - log(-t) + / log 


ICO 


A^ - cr" 
-t 


dv{a) , 


where (j = 0,..., s — 1) are the components of C (as defined in the proof of Lemma 1.3). Since on each 
C-l we have that dO is normalized to - we obtain 


lim — logp„(A) = [ log(A® — cr^)(ii>(cr) = / log(A — crwy(ii>(cr) 

n-i-oo n Jco ^ Jco 

= V / log(A - a)di)(a) = [ log(A - a)dv{a) , 

.■_n Jci Jc 


with u> = and where in the last steps we used the symmetry of dv. Hence we have obtained the 
relation 

lim -logpn(A)= / log(A - cr)df'(cr), 

n—>-oo n 

uniformly for A in compact subsets of unbounded component of C\C. Furthermore C is the boundary of 
its polynomial convex hull. Then it follows that the measure v is the weak-star limit of the zero counting 
measure Vn of the polynomials Pn(A) ioi n = kd 1, I = Q,s — 2 (see [44] and [49] Chapter 3). The 
proof of Theorem |1.4| is then completed. 
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